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Abstract 

The relevance of the dipole configurations of quarks in forming nucleus structure func- 
tions is discussed. It is shown that a radiation generated by dipole configurations while 
moving relativistically along their axises is described by distributions which are finite and 
infrared stable in low transverse momentum region. It is argued that there is an expo- 
nential transition to the perturbative regime of large transverse momenta and its power 
is defined by the distance between the dipole charges in its rest frame m^ 1 . 



1 Introduction 



Theory of dense relativistic gluon (parton) systems is recently recognized to be crucially 
important solving many (even long time pending) problems of high energy hadron physics 
. But on eve of new round of the experiments with ultrarelativistic nuclear beams (at 
RHIC and LHC) this subject becomes a real challenge because of an extreme necessity 
to answer most reliably the question about the initial conditions and early evolution 
of the quark-gluon plasma (if it is formed) as emphasized in |^]. Apparently clear, an 
answer depends essentially on the behaviour of the nuclear structure functions which are 
calculable perturbatively in QCD Q. However, extrapolating these results even to the 
region of moderate parton transverse momenta takes very special efforts H and leads to 
the singular distributions for the small momentum values. 

Moreover, a rather provocative idea in this context launched in []|] suggests for the 
largest nuclei at very high energy the initial parton density could become so high that 
the intrinsic transverse parton momenta reach the magnitudes where an overlap with 
the region mastered by perturbative QCD already occurs. It treats the structure func- 
tion of nucleus as an aggregate of quantum fluctuations on the ground of classical gluon 
fields generated by ultrarelativistic (anti-)quarks and argues that respective classical beam 
bremsstrahlung could modify many results which are discussed as possible quark-gluon 
plasma signatures. In the course of such calculations 0, || light cone QCD is handled 
in the light cone gauge @ because of taking the limit of v — > 1 becomes quite nontrivial 
problem || and the corresponding potentials turn out pretty singular and having no the 
direct physical meaning. 

In the present paper dealing with classical electrodynamics in order to make the re- 
sults somewhat easily understandable we show that the field strengths of particular dipole 
configurations regularized by averaging E reg , H reg justify taking the immediate limit of 
v — * 1 for physical observables and provide reasonable asymptotic behaviour for the struc- 
ture functions. We believe it could be indicative to appraise nonperturbative region role 
contributing to the behaviour of nuclear structure functions in the McLerran-Venugopalan 
approach. 

The paper is simply structured. In Sect. 1 we mainly remind the basic formulae for 
the potentials and field strengths of ultrarelativistic particles in classical electrodynamics. 
The radiation of a fast-moving classical charge in the fields considered is calculated in 
Sect. 2. We conclude discussing the relevance of these results to describe the nuclear 
structure functions in the final Sect. 3. 

2 Field of ultrarelativistic particles in electrody- 
namics 

If one considers the field of a classical charge e freely moving with the constant velocity v 
the following expression for the potential is valid demonstrating obviously the cylindrical 
symmetry of a solution in the Minkowski space (c = 1) 

M*) = ^ (!) 
K 

where u u = (uo, u), u = vuqti, uq = ^„ -. /9 and n is the unit vector along particle 

(1 — v ) 1 

velocity direction whereas R 2 = x 2 + (xu) 2 = x\ + ((xn) — vxq) 2 /(l — v 2 ). Then the field 
strengths are given by 
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jg(.r) = 6 y° U \ H = vxE. (2) 



At v = 1 the potential does exist but is singular in the plane R = 0, in particular, 
moving along z-axis it happens to be that 
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\Z - X \ 

A X = Ay = 0, 

A,= ' 



- xq\ 

The field strengths disappear (equal to zero) in the limit v — > 1, besides the plane 
z = xo where the transverse components tend to blow up 

x y 
E z = 0, E x = e uq <r ■ E y = e u 

here with p 2 = x 2 + y 2 . Actually, nonexisting finite limit of physical observable quantity, 
like the field strengths in this case, signifies simply that a real physical charge can't be 
accelerated up to the light velocity because of infinite amount of energy necessary to do 
so. 

The potential of randomly distributed charges with the center of mass ultrarelativisti- 
cally moving in the direction of vector n is given by the following integral over the charge 
density distribution in the rest frame of q(z),i.e. 

A^x) = J dz q(z) A^(x,z), 

An(x, z) = , M - , 
R (x, z) 

where R 2 = — Xq + (x — z') 2 + [—xq uq + (x — z' ,u)] 2 , and the vector z absorbs the 
Lorentz contraction of the system of charges in the factor Uq 1 in the motion direction, 
thus z = (zn)n/uQ + z — (zn)n. 

In the particular case of a dipole with charges e and — e being placed at the points 
z\ = —d/2, Z2 = d/2 of the dipole rest frame, the electric field components become 



E± = e Uq 



x ± + d ± /2 x ± -d ± /2] 

Ri Ri J 

(dn) 1 e f (dn) 



e ( (dn) 1 e f 

II = 773 1 u o( xn ) ~ vu x H — > ~ —r^i u (xn) - vu x 

Ri I 2 ) /£ 2 I 



(4) 



with R\ = R'(zi), z = l,2. 

Although these fields are singular in the limit v — ► 1, nevertheless one may notice their 
orthogonal (relative to the motion direction) components are developing two symmetric 
peaks opposite-orientated at = . Let us trace now the reaction of fast moving 
particle passing through the field of such a configuration. Due to the Lorentz contraction 
this particle is exposed to two very short <5-like pulses of opposite signs. If the contraction 
is large enough that the time interval between two pulses rj ~ T is less than the classical 
radius of particle a = e 2 /m, d ~ a with mass m, there is no enough time to respond both 
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separate pulses and the charge is sensitive to the smeared (averaged in time) field only. 
Indeed, the evolution of charge radiation is controlled by the Newtonian equation with 
the radiative friction 



w = g(t) + 2e 2 w/3m, 
and its general retarded solution is 



w = WQ 



+ f g{r)dT+ f%(r)e(*- T 

J — OO J t 



if g(t) is the impulse generated by an external field, k = 2e 2 /3m. When two <5— like contra- 
peaks of the amplitude ~ G separated in time ~ T at the condition T/k <C 1 are present, 
the particle velocity w differs from the initial value wq by the quantity of G T 2 / k 2 order 
(in the order of T/k the pulses cancel each other). Since G ~ vuqGo and T ~ Tq/vuq 
with Go and To being fixed in the rest frame, then in the limit v — * 1 it occurs that 
G T 2 1 k 2 — > 0. Hence, it is well grounded to write down for the general solution 



w = w + (g(r))dr+ {g(T))e^/ K dr, 



where (g(r)) is an impulse smeared in time. 

Thus, as the regularized field we should treat one averaged over the interval [xm, xof] = 
[{xn)/v — (dn)/2vuQ, (xn)/v + (dn)/2vuo], i.e. defined by 

w = s r e dx °- ^ 

{dn) J X0 , 

In particular, if the charge distribution takes such a form that (d±) ~ (averaging over 
charge density), then we have 

^ll) = 7^lurT- ua *\ 2U/2 j. (E±)*Q. (6) 
(dn) y\x±\ [{dnf + xj} 1 -'^ J 

It is interesting to notice that our procedure of regularizing by averaging is similar, 
in a sense, to the potential regularization as done in Ref.||] dealing with the light-cone 
variables and then widely used ||,||,[|l(J. If we integrate now the transverse field of a 
singled charge over dz (when moving in z-direction) on the infinitesimal interval [zi, zA = 
[z — vxq — b/ua, z — vxq + b/uo] (v — > 1) 

'f , e x±_ 2b 
E ± dz 



p 2 ( 6 2 +p2) l/2- 

at large enough value of b 3> p then it is valid 

lim f E±dz — > — tt^, 
p/b->Q J Zi - p 2 

i.e. the transverse field may be presented in the form of 5 - function E± = ^ e ^' L S(z — xq), 

what is equivalent to the equation (8) of Ref. Q. It was calculated there starting from 
the initial potential of Eq.(l), omitting the retarded part (x + = (z + xq)/2 ~ 0) and 
performing the singular gauge transformation of the longitudinal potential components to 
the transverse ones. Strictly speaking it's not permissible to drop out the x+-dependence 
because the corresponding terms are of the same order in (1- v) as ones used in the further 
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calculations of Ref. 0] (the entire transverse field is zero in this approach as well if we 
consider the dipole configurations of longitudinal orientation). For the transverse dipoles 

^ o x( J x ± + d ± /2 x ± -dj_/2 \ 
E ^ = 2e5{Z -^{ { X , + d j2 ? - {x ,- d j2)A 

and we obtain the entire transverse field equal to zero again but after averaging over its 
orientation in (x,y)-plane for a homogeneous (like in KR of Ref. §) distribution. 

3 Radiation of fast charged particle in averaged 
field of relativistic dipole 

The radiation energy of a fast moving classical charge in an external electromagnetic field 



is determined as [11] 



3 m 2 1 — vi 2 

where e\ and V\ are the charge and velocity of radiating particle, respectively. Here we are 
interested in the situation when a particle and an ultrarelativistic dipole oriented along its 
motion direction are going towards each other. In accordance with previous consideration 
we should substitute the average field instead of instant values of field strengths in the 
radiation formula. Such replacement results in disappearing the singular part of radiation 
from the transverse components of field, because (E±) ~ 0, (H±) ~ 0. Then the same 
factors (1 — v 2 ) available in the numerator and denominator are cancelled out and the 
energy radiated per unit time comes about to be independent of the particle energy 

dAe _ 2e\ 2 

~dT - 3^ W • (7) 

here we neglect, of course, a particle influence on a dipole in the first approximation. 

For comparison we give the radiation intensity of a fast charged particle in the field of 
immobile dipole with its orientation along the direction of particle propagation in eikonal 
approximation: 

3 vr ef e 2 d 2 , n 5 , 

Ae = i— , — 3 + * ) 

64 m z p b v 1 1 — vf 

(see also a similar exercise for the radiation in the field of Coulomb center [|ll|] ). This 
radiation is composed by two different components, i.e. the singular part of transverse 
components and the regular one of longitudinal components. In fact, we can neglect the 
longitudinal radiation in the limit v\ — > 1. Our averaging procedure and substitution of 
(E), {H) in radiation formula instead of E and H kill a singular radiation and we must 
take into consideration small longitudinal radiation 'surviving'. This substitution arises 
as a leading term of perturbative expansion in T/k which we are able to argue basing on 
the Newtonian equation (though it can be done in the framework of relativistic dynamics) . 
Any other terms should be omitted in the solution of dynamical equation with radiation 
force included. Let's trace now the consequences of averaging procedure. 
Integrating Eq.(7) over the impact parameter we have 



dAe _ f d 2 k L 2ef 
~dT ~ J (2vr) 2 3m 



2 EfAk±), (8) 
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where 



E ]l (k±)= f dx ± e~ i{k±x±) (E\\ 



Aire 
d k i 



(9) 



and Ki/ 2 (z) is the modified Bessel- function. Then asymptotic behaviours of the field are 
the following 

lim Eu(k ± ) -> 47re , lim Siiffcj.) -» 4t^{1 - e" d fc± i . 
fc x ^o 11 fex^oo 11 a k± 

It results in the finite value of the distribution of radiated photons as obtained from 
Eq.(8) divided by the photon energy ko, and being proportional to the Fourier component 
of electric field squared, of course, in the limit k± — > . It approaches exponentially 
the perturbative behaviour ~ \jk\ at the transverse momenta while large enough. Let's 
emphasize the rate of transition to the asymptotic regime is regulated by the distance 
between the dipole charges in its rest frame ~ to" 1 and the radiation energy is proportional 
to (e 6 ). 

The x distribution (x is a portion of longitudinal momentum carrying by an individual 
parton) takes well-known 'soft' form dk^/ko « dx/x, because 8(z — xq) ko = k% k\. 
In order to get beyond this approximation we need to keep the next order terms of the 
(1 — w 2 )-expansion. Moreover, in lieu of the coherent summation of the fields as in Eq.(5) 
one should average them taking into account the phase shifts between different points. 
Calculating with the corresponding phase factor taken (moving along the z-direction) 
leads to 

fi + 4) 1/2 -i-^ 

{E x e"») = 6 - 2ik 3 e*™ x ± ^ [/ + ^ ' , 



(E z e ik3Z ) = 2 e ik3VX ° 



1 1 



P [p' + d 2 ] 1 / 2 



These expressions make it clear that the transverse field (unlike the longitudinal one) may 
be dropped out when one estimates the radiation behaviour in low transverse momentum 
region (large impact parameters p). Although the transverse field doesn't identically 
disappear and allows us to take the limit v — * 1, it's justified, nevertheless, because of the 
additional suppression d 2 j p 2 present comparing to the longitudinal field. The coherent 
regime is associated with the constraint k%d < uq , which is apparently valid in the limit 
of v — > 1 for any k%. We understand, however, these delicate problems require the full 
detail investigations in future. 

In fact, the analysis performed may be extended to non-abelian theory. As shown in 
Ref. |l2j at the small enough coupling constant g 2 /An < 1 non-abelian theory is steady 
similar to the electrodynamics as to the problems of both dipole and induced radiation 
description. The influence of non-linear terms, in particular, of 'three gluon vertex' mainly 
discussed, is reduced for the dipole configuration to provide only pure rotation of the 
vectors of particle colour charges in the 'isotopic' space (for SU(2) non-abelian group) 
which is certainly unobservable and doesn't lead to any physical effects. What one should 
modify only is to replace the charge product e\ e 2 of the interacting particles by the 
scalar product (PiP 2 ) of the vectors of particle charges in the 'isotopic' space and to take 
also into consideration the energy of gluomagnetic field generated by the chromoelectric 
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charges Pi and P2 ■ Thus, the existence of additional distribution over the relative angles 
of 'isotopic' vector charges of quark-antiquark pairs in nuclei foresees major distinction of 
non-abelian theory from electrodynamics in such an analysis. 



4 Conclusion 

Using the classical electrodynamics as an example we have demonstrated above the limit of 
v — ► 1 does exist for the physical observable quantities, what are the average field strengths 
of dipole configurations oriented along their propagation direction ('longitudinal' dipole), 
and this qualitative picture does not change when non-abelian theory is considered. In 
a sense it justifies some intuitive representations about the structure of classical gluon 
fields generated by ultrarelativistic nuclei collision and about nuclear structure functions 
developed by McLerran and Venugopalan in ||. The gluon transverse momentum dis- 
tribution here obtained agrees qualitatively with the distribution presented in Fig.l of 



Ref. [1C] though essential difference is that we obtain the infrared stable and finite result 
at k± — > and the transition to the perturbative regime occurs at the transverse momenta 
k\ ~ d~ 2 ~ m 2 . We are free to take the minimal size of hadron as the mean distance 
between dipole charges in its rest frame. Unfortunately, we are unable now to point out 
the well grounded mechanism to orientate the colour dipoles along their propagation but 
we believe (and may argue) the rather plausible hypothesis is that the dipoles are aligned 
along the external accelerating field because such configurations in electrodynamics are 
energetically more favourable. 

Let us summarize now our main results. We suggest the regularization procedure of 
classical fields as being generated in ultrarelativistic heavy ion collisions and relate the 
field smearing to the Lorentz contraction of field configurations. If characteristic field 
configuration 'size' is less than the classical particle radius then this particle is able to 
'feel' an average field only. Operating with the physical observables we have immediately 
calculated the average field of fast moving 'longitudinal' dipole (Eq. (6)) in the limit v — ► 1 
and then the radiation energy of a fast moving charge in this field (Eqs. (8), (9)). The 
transverse momentum distribution of radiated photons happens to be finite at small values 
and exponentially runs into the perturbative regime ~ l//cj_ at larger k\. Applying proper 
quark-antiquark configurations for the calculation of structure functions leads to the same 
result as for the photons distribution and agrees qualitatively with result obtained earlier 



in Ref. [10|. Unlike the light-cone method every stage of our approach is physically 
meaningful because we are dealing with observable values only. Seems, including next 
order terms of expansion in (1 — v 2 ) could illuminate the origin of structure function 
behaviour at small x what we see as the subject of our further investigation. 

The work has been initiated by numerous discussions of the McLerran- Venugopalan 
model with V.Goloviznin, M.Gyulassy, Yu.Kovchegov, A.Kovner, A.Leonidov, A.Makhlin, 
L. McLerran, D.Rischke, H.Satz and R. Venugopalan. The financial support of RFFI 
(Grants 96-02-16303, 96-02-00088G, 97-02-17491) and INTAS (Grant 96-678) is greatly 
acknowledged. 
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